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A Green functions approach is used to study superconductivity in nanonhns and nanowires. We 
show that the superconducting condensate results from the multimodal entanglement, or internal 
Josephson coupling, of the subcondensates associated with the manifold of Fermi surface subparts 
resulting from size-quantisation. This entanglement is of critical importance in these systems, since 
without it superconductivity would be extremely weak, if not completely negligible. Further, the 
multimodal character of the condensate generally results in multigap superconductivity, with great 
quantitative consequence for the value of the critical parameters. Our approach suggests that these 
are universal characteristics of confined superconductors. 

PACS numbers: 74.20.-z, 74.78.-w, 74.81.-g 



In recent years, superconductivity has been stud- 
ied with great interest in nanoscale systems, such as 
nanofilms (NFs) |l|-l4( and nanowires (NWs) [5|, [6[ . From 
the theory point of view, there was interest in su- 
perconductivity in NFs [7j long before the advent of 
nanoscience. The great advances in materials synthesis 
technology, however, allows today the fabrication of high 
quality nanoscale samples, in which the effects of confine- 
ment, or quantum size effects, can be closely examined. 
Thus, for instance, the oscillations of the superconduct- 
ing gap with film thickness predicted in Ref. 7 have now 
been observed convincingly in experiment, albeit with a 
weaker amplitude than the one predicted Recent 
theoretical work has shown similar quantum size-induced 
oscillations in the critical temperature, specific heat, and 
critical field in NFs |1, Q as well as in NWs [To|, EL 
Interestingly, it is found that, for sufficiently small con- 
finement length (i.e., film thickness or wire radius), the 
oscillations can drive the critical temperature well above 
the bulk value. This was not observed in the ultrathin 
NFs of Refs. [ij-y, but appears to be the case in the NFs 
of Ref. [H, and also in NWs @, @. 



Here we show that confinement has far reaching conse- 
quences regarding the character of the superconducting 
condensate itself. In this regard, the study of supercon- 
ductivity in nanosystems is of broader interest. Indeed, 
superfluidity also fascinates researchers in atomic physics 
[13| and nuclear physics [14j|, where confinement occurs 
only naturally. Superfluidity is a macroscopic quantum 
phenomenon, in which the order parameter, or conden- 
sate "wave function" plays a central role. Thus, much 
effort is put in trying to characterise it in the different 
systems in which superfluidity is observed [l3|, 15]. In 



this work, we use a Green functions approach to super- 
conductivity in nanosystems, bringing to light previously 
unrecognised but significant confinement effects. We find 
that the splitting of the Fermi surface in a discrete set 
of nonintersecting parts because of the size-quantisation 
of the energy levels results inevitably in a condensate of 



composite nature. That is, the condensate arises from the 
multimodal entanglement, or internal Josephson coupling 
[la LL7j , of the subcondensates associated with the Fermi 
surface subparts. Further, the entanglement is crucial, 
because without it the systems collapses into uncorre- 
cted subcondensates, in which superconductivity is dra- 
matically weaker, if not completely negligible. Finally, 
the manifold of subcondensates will generally result in 
multigap superconductivity, with significant impact on 
the predicted value of the critical temperature, compared 
to single-gap models. The generality of our formalism 
suggests that these characteristics are universal proper- 
ties of confined superconductors. 

Consider a system of quasiparticles with a weak 
attractive effective interaction coupling only particles 
with opposite spin (e.g., the net effect of phonon ex- 
change and the screened Coulomb interaction), which 
will eventually couple only particles near the Fermi 
level (chemical potential), \i. We have the Hamil- 
tonian H = Ho + Hi — uN : where Ho describes 
the noninteracting system, N is the number operator, 

and Hi = i/dVdV^IW^^rOveff^rO^ICrO^tW- 
Consider next the Green function, G(yy' ,t) = 

-(T r ^{(rr)^ t (r / 0)). As in BCS theory [18], to solve the 
equation of motion for Q we introduce the Gorkov func- 
tions F(ty',t) = -(T r ?/i t (rr)^(r / 0)) and P (rr',T) = 
— (T r ?/i|(rr)^(r / 0)), and take the mean- field ap- 
proximation (T T $1 (r i n ) ^ (r 2 r 2 ) ^ (r 3 r 3 ) ip j (r 4 r 4 ) ) ~ 
— J r (r3r2,T3 — r2)J r ^"(rir4,ri — T4). The resulting cou- 
pled equations for Q and read, in frequency domain, 

£ p (r)£(rr',u; p ) - j d 3 zA(ra)^(xr',a; p ) = H5(r - r'), 
C p {y)P{yy' ,u p ) - J d 3 £A(rx)£(xr', uj p ) = 0. (1) 

Here, C p (y) = ihu p — Hq(y) + /i, with uo p is a fermionic 
frequency, denoting a r-Fourier transformed function, 
and we have introduced A (it") = v G s (r, r /, )J 7 (rr ,/ , 0), 



2 



i.e., a nonlocal pairing potential [14 Il9|. 

Assume now that Hq is such that its eigenstates form 
an orthonormal set, i.e., H$\v) = E v \v), with (v\v') — 
5 UyU > . Inserting the expansions i/j a (rr) = ^2 U ^ u (r)c U(T (r) 

and i>l(rr) = Z^^( r ) c L( T ) in Ec L s - O leads to [ witn > 
e.g., £(rr> p ) = Mr)W>P)GW,u p )] 

[zfi^p - e^] G{vv' ,uo p ) + ^ A(^ft)J^~(ftz/,cap) = fkW, 
+ ^A*(i/k){7(ki/ / ,w p ) = 0. (2) 

ft 

Here, e v = E v - n and A(W) = V^^/J 7 ^', 0), 
with Vw' ,kk' == — (^^'l^eff |^^')- The f° rma l solution of 
Eqs. ([2]) can readily be written in closed form. But we 
further assume that v G s couples only time-reversed states 
[20|. If | — v) and \v) denote time-reversed states, we 

haVe VvV f ,K,K f ~ ^V—V^K—K^—VV'^—KK' — V VK b— v v l 6— K K ' 

and A(W) = 5- vu ' V^J-^— ft 0) = A(v)5-„„>. 
As in BCS theory, it is straightforward to deduce 
that G(vv',u p ) = -5 vv >{iuj p + e„/h)/(u 2 + Q/h 2 ) and 

P(vv'^ v ) = 5^,(A*(is)/h)/(u 2 p + with H = 

e 2 + |A(V)| 2 . Thus, the gap values are given by 
the coefficients of the expansion of A(rr // ) over the 
quasiparticle states. Finally, noting that J r (— kk,0) = 
F(—kk,u p )/ f3h, we obtain the gap equation 



A(z/) = ^K^A(z/ , )^tanh 



2fc B T" 



(3) 



Once A(^) is determined, the condensate wave function, 
defined by ^/(rr') = ^"(rr',0) [21], and the pairing po- 
tential can be calculated from 

*(rr') = ^V,(r)V-,(r')AM^-tanh^, (4) 
A(rr / )=^^WV'-,(r')AM. (5) 

In the following we apply our formalism to NFs and 
NWs. We take parameters corresponding to Al (22| . 
which is a weak coupling superconductor (so a mean- 
field approach is applicable), and is also free-electron-like. 
Given that the spatial dependence of v G s is not really 
known, it is best to use phenomenologically motivated 
V^/'s in Eq. (j3j) (thereby, implicitly defining the spatial 
form of f e ff). We use as reference the BCS coupling con- 
stant for the bulk material, Vb, estimated from the ex- 
perimental value of T c and k B T c = 1.13/kJ D e- 1/iVoVo pjjj. 
Below, unless otherwise stated, energies are in Rydbergs 
(Ry), and lengths in a$. 

To model a NF we follow Ref.[24|. Briefly, the system of 
quasiparticles is in a potential well defined by two large 
planes of side L, a distance d apart, with V(p, z) = 
for < z < d, and oo otherwise. The quasiparti- 
cles states are -0 qn (p, z) = (2/L 2 d) 1 ^ 2 e' lCL ' p sin a n z, where 



a n = nn/d and q = 27r(/,ra)/L, with Z,ra G Z. Thus, 
1 1/) = |qn) and | — v) = | — qn) are time-reversed states. 
The energy levels are given by £? qn = q 2 + a 2 (for sim- 
plicity, the quasiparticle mass is taken equal to the bare 
electron mass). The Fermi "surface" breaks into a set of 
concentric circumferences of radii = (/i — a 2 ) 1 / 2 . For 
the V vv > in Eq. (|3]) we take a BCS-type model [l8|, 



qn,q'n' 



L 2 d 



#(e w - |e qn |)#(e w - |e q / n /|), (6) 



e w defining the energy window around \i within which 
v e ff is effective. We estimate the U n71 > with a contact 
potential v e ^(r — r f ) = — Vb5(r — r') [l8|,|25[, hence U nn > — 
Vb(l + <W/ 2 ) [24]. Note that Eq. © will lead in general 
to a multigap equation, similar to the expression of Suhl 
et al. [26(| . But because the contact interaction results 
in off-diagonal U nn ' 's that are all equal, there is only one 
gap. Thereafter, it is straightforward to derive the results 
of previous authors (cf., e.g., Refs. 0, 0). 

It is important to recognise, however, is that even with 
the simple contact potential, this still is a system with 
multiple subcondensates. To see this, we rewrite the 
Hamiltonian in second-quantised form. In close analogy 
to the two band case studied by Leggett [16j, one finds 
H = H D +Hj, with 
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n qcr 



U n ^ C qnt C -qn| C -q / n| C q / nt 



qq' 



^ J — [ J ^ C qnt C -qn4, C -q / n / 4 r C q / n / t * ^ 



qq' 



Here, is the Hamiltonian of the independent conden- 
sates, with U n = 3Vo/2, while Hj represents an internal 
Josephson coupling [l6|, with J = Vq. Thus, the conden- 
sate is given by a multimodal entanglement of subcon- 
densates [17]. The entanglement is beautifully illustrated 
by the resulting interference pattern in the probability 
density, ^(rr')! 2 . In the present case, Eq. (|4]) leads to 
* = *(p, z- z') = En *n(p, Z] z'), with (at T = K) 



4ir 2 d 



I n (p) sin a n z sin a n ^, (8) 



where I n (p) = f' dqqJo(q\p\)/£q n [27]. Given two quasi- 
particles of opposite spin, at (0,z f ) and (p,z), respec- 
tively, \I/(p, z\z') is their pairing probability amplitude. 
In Fig. Ufa), we plot | ^ n * n (p, ^)| 2 for (p, 2) in the 
xOt/ plane, for a d — 10 ao NF, and contrast it to 
Y, n \^n{p^-z')\ 2 [cf. Fig. 1Kb)], to highlight the inter- 
ference effects. To choose z' \ we calculated first the local 
pair density, g s {z') = J d 2 p dz\^(p, z; z')\ 2 . One readily 
finds Qs(z') = (A/167r 3 <i) arctan(e w /A) sin 2 a n 2; / . In 
Fig.E^a) we plot g s (z') for three d values. For d = 10 ao, 
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|^(z'=0.79d)| 2 (a Q 6 ) -13 Z n |^ n (z'=0.79d)| 2 (a Q 6 ) -n 
xlO n n °xlO 
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FIG. 1: (color online) Condensate probability density (a) 



|£„<MP.*;*')r and £J*„(P, 



(b), plotted in the 



xOz plane, at z — 0.79d, for d — 10 ao (see text). In (b) 
there are no interference effects. 




0.5 0.75 
z7d 



0.25 0.5 0.75 
Factor f 



FIG. 2: (color online) (a) Local pair density, g s (z'), for three d 
values (in units of the bulk electron density, n). The number 
of maxima indicate the number of subcondensates. (b) A 
and T c change strongly with factor / (see text). At / = 
(decoupled condensates), T c ~ 10" 3 T c b , and A ~ 10" 3 A b . 
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to a discrete set { — k^ nn \k^ nn) } mn . The energy bands 



l /R , and the Fermi surface reduces 

are now 1-dimensional, while they were 2-dimensional 
in the NFs. This gives rise to important quantitative 
differences between the two cases regarding the be- 
haviour of their properties as a function of confining 
length [HI [H|. Here we focus, however, on the multigap 
character of superconductivity in NWs. To see this, let 
us approximate the V vv > in Eq. (j3]) by 



Qs(z') is maximum at z' = 0.79a 7 (note that the number 
of maxima in q s (z') corresponds the number of subcon- 
densates in the film). 

Furthermore, the strength of the J coupling is of crit- 
ical importance, its magnitude largely determining the 
value of the critical parameters. Indeed, for a renor- 
malised coupling J = /Vb, with / < 1, both A and 
T c fall dramatically as / decreases. We illustrate this in 
Fig. Efb), for d = 10 ao- At / = 1, the critical parame- 
ters are significantly higher than the bulk values, namely 
T c jT h c = 2.51 and A/A 6 = 2.60. For / = 0, i.e., decou- 
pled condensates, A/ A 6 and T c /T c 6 are negligible, of the 
order of 10 -3 . In contrast, in Refs. [H-E!, A and T c are 
found to be a large fraction of the bulk values, requiring 
0-5 < / < 1 in our model, i.e., a substantial coupling. 
A value / < 1 is easily understood, since interband pair 
scattering requires a minimum momentum transfer, so 
has a smaller scattering phase space volume than intra- 
band scattering [an aspect not accounted for in Eq. flgj)]. 
In fact, this may be another reason why in experiment 
the critical parameters are lower than in the bulk. 

We now turn our attention to NWs. The quasiparti- 
cles are now in a cylindrical potential well of radius R 
and length L: V(p, </>, z) =0 for p < R, and oo other- 
wise. The quasiparticle states are [To| ipkmnip, <!>•> z) = 
M^^l+i^)]" 17 ^^!^^/^)^^^, where 
J m is the ra-th order Bessel function of the 1st kind 
and r] mn is its n-th zero [lij, and k — 2nl/L, with 
I G Z. Here, \v) = \kmn) and | —v) = | —k —rn n) are 
time-reversed states. The eigenenergies are given by 



Vkmn,k'm'n' — o r-» r @(^w \^-kmn\)^{^-w \^k'm'n' |)- 
7T Z RL 

(9) 

To estimate the /7 mn5m / n / we again use a contact poten- 
tial. Unlike the NF case, the off-diagonal elements are 
different from each other. This immediately results in 
multiple gaps, A mn . In Fig. [3fa) we plot the gap values 
as a function of temperature for a R = 7.5 ao NW [29| . 
In this case there are seven occupied bands, thus seven 
subcondensates. The A mn (0) values depend on the inter- 
play between the Z7 mn5m / n / strengths and how far from p 
are the bottoms of the bands [recall that in 1-dimension 
the density of states has an (integrable) singularity at 
k = 0]. T c and the A mn oscillate strongly as a function 
i?, rising sharply when the bottom of a newly occupied 
band falls below p as R increases. This is illustrated for 
T c in Fig. [3fb) (upper panel) [30]. Although similar to 
the oscillations found in single-gap models [To|, QjJ , the 
multigap character of the condensate results in signifi- 
cant quantitative differences. Indeed, Fig. [3^b), lower 
panel, shows the plot of the ratio of T c 's obtained in the 
multigap and single-gap cases (the latter, T c s ^, is obtained 
by approximating the U mnjrn f n f by their average value, 
Umn,m f ri)- We see that T^ 9 can be more than 100% too 
low respect to the multigap value. As one would expect, 
the magnitude of the J coupling is just as critical here 
as in NFs. Indeed, setting /7 mn ^ m / n / = in the gap 
equation results in uncorrelated condensates, with x^ mn ^ 
values largely reduced respect to the true T c . For exam- 
ple, in a R = 5 a NW, for which T c = 2.12 T c 6 , there 
would be three condensates, with critical temperatures 
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FIG. 3: (color online) (a) Plot of the seven A mn (T) in a 
R = 7.5 do NW; T c — 4.2 T c b . (b) Upper panel: T c as a function 
of R. T c increases sharply when a new band starts to be 
occupied. The horizontal line indicates T c b . For large R, T c 
tends to T c b (not shown here). Lower panel: The ratio of T c to 
the single-gap value, T c s ^, shows that they differ significantly, 
(c) The pairing potential, A, and order parameter, ^, are not 
proportional to each other (here shown in the r = r' limit). 



T c (21) ~ 0.078 Tj\ T c (11) 



10- 5 Tj\ andT c 
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We add that, because the matrix elements decrease 
with confining length, a finite J coupling is essential to 
obtain the bulk values of the critical parameters in the 
limit of large systems (i.e., R — > oo in NWs and d —> oo 
in NFs). Also, as defined in Eqs. (j4]) and (j5j), the pairing 
potential and the order parameter are not proportional 
to each other (unlike in homogeneous systems [18]). even 
in the r = r' limit. For example, in Fig.[3fc) we compare 
A(p) (renormalised, for comparison) and ^f(p) (in that 
limit both depend only on p) for the R — 7.5 ao wire. 
So our A(p) is not equivalent to the "order parameter" 
in other approaches [To|, [H|. Also, our A(p) should not 
be confused with the spatially varying gap seen, e.g., in 
some high-T c superconductors [31]. Indeed, in our case 
the gap(s) are constant throughout the system. 
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